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Abstract: In the first part of this paper, we outline the construction of an inflationary 
cosmology in the framework where inflation is described by a universally evolving scalar 
field (p with potential V{(j))- By considering a generic situation that inflaton attains a 
nearly constant velocity, during inflation, mp^\d(j)/dN\ = a + /3exp(/3A^) (where = Ina 
is the e-folding time), we reconstruct a scalar potential and find the conditions that have 
to satisfied by the (reconstructed) potential to be consistent with the WMAP inflationary 
data. The consistency of our model with WMAP result (such as = 0.95llQQjg and 
r < 0.3) would require 0.16 < a < 0.26 and (3 < 0. The running of spectral index, 
a = dn^/dln k, is found to be small for a wide range of a. 

In the second part of this paper, we introduce a novel approach of constructing dark 
energy within the context of the standard scalar-tensor theory. The assumption that a 
scalar field might roll with a nearly constant velocity, during inflation, can also be applied 
to quintessence or dark energy models. For the minimally coupled quintessence, Oq = 
dA{Q)/d[KQ) = (where A{Q) is the standard matter-quintessence coupling), the dark 
energy equation of state in the range — 1 < Wj^^ < —0.82 can be obtained for < a < 0.63. 
For a < 0.1, the model allows for only modest evolution of dark energy density with 
redshift. We also show, under certain conditions, that the Og > solution decreases the 
dark energy equation of state Wq with decreasing redshift as compared to the = 
solution. This effect can be opposite in the Oq < case. The effect of the matter- 
quintessence coupling can be significant only if la^j ^0.1, while a small coupling |ag| < 0.1 
will have almost no effect on cosmological parameters, including Qq, Wq and H{z). The 
best fit value of Uq in our model is found to be Uq ~ 0.06, but it may contain significant 
numerical errors, viz Uq = 0.06 it 0.35, which thereby implies the consistency of our model 
with general relativity (for which Og = 0) at Icr level. 
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1. 


Introduction and Overview 




It 


is true and remarkable that our understanding 


; of the physical universe has deepened 



profoundly in the last few decades through thoughts, experiments and observations. Along 
with significant advancements in observational cosmology [1-3], Einstein's general relativity 
has been established as a successful classical theory of gravitational interactions, from scales 
of millimeters through to kiloparsecs (1 pc= 3.27 light years). It has also been learned 
that at very short distance scales large quantum fluctuations make gravity very strongly 
interacting, implying that general relativity cannot be used to probe spacetime (geometry) 
for distances close to Planck's length, Ip ~ 10~^^ cm. In addition to this difficulty, three 
striking facts about nature's clues suggest that we are missing a few important parts of the 
picture, notably the extreme weakness of gravity relative to the other forces, the huge size 
and flatness of the observable universe, and the late time cosmic acceleration. 

Much is not understood: what is the nature of the mysterious smooth dark energy 
and the clumped non-baryonic dark-matter, which respectively form 73% and 22% of the 
mass-energy in the universe. That means, we do not see and really understand yet about 
95% of the total matter density of the universe. To understand the need for dark energy, or 
a mysterious force propelling the universe, and dark matter, one has to look at the different 
constituents of the universe, their properties and observational evidences (for reviews, see, 
e.g. [4-7]). The current standard model of cosmology somehow combines the original hot 
big bang model and the early universe inflation, by virtue of the existence of a fundamental 
scalar field, called inflaton. The standard model of cosmology is, however, not completely 
satisfactory and it appears to have some gaps. If the universe is currently accelerating (on 
largest scales), what recent observations seem to indicate, then we need in the fabric of 
the cosmos a self-repulsive dark energy component, or a cosmological constant term, which 
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had almost no role in the early universe, or need to modify Einstein's theory of gravity on 
largest scales in order to explain this acceleration. 

When in 1917 Einstein proposed the field equations for general relativity 

RfMi^ - \ Rgfiu + A^Mi^ = SttGjv Tf,^, (1.1) 

he had the choice of adding an extra term proportional to the metric g^iy either on the left- 
hand or right-hand side of eq. (1.1). This extra term, so-called the cosmological constant 
A, is not fixed by the structure of the theory. One also finds no good reason to set it to 
zero either, unless the underlying theory is purely super symmetric. Adding the term Ag^i, 
on the left hand side of his famous equation, Einstein used to tune the constant A in such a 
way that he would get a non-expanding solution. Einstein later dismissed the cosmological 
constant as his "greatest blunder", when Hubble found a clear indication for an (ever) 
expanding universe. Today this constant is mainly written on the right hand side of the 
Einstein equations but still with a positive sign, which therefore acts as an extra repulsive 
force (or dark energy) in cosmological (time-dependent) backgrounds. 

Before presenting further thoughts on the nature of this puzzling form of energy, it 
is logical to recapitulate the independent pieces of evidence for its existence. The key 
measurements, leading to the result of DE density fraction being f^DE — 0.7 have been 
made, rather unexpectedly, in 1998 by two independent groups (Supernova Cosmology 
Project and High-z Supernova Search Team) [1]. These observations revealed, for the first 
time, that the universe is not only expanding now but its expansion is speeding up for the 
last 5 — 6 billion years, i.e. since when the redshift z dropped below 0.85. 

Evidence for the existence of dark energy also comes from observations of the Cos- 
mic Microwave Background (CMB) for which the most recent ones have been obtained 
by NASA's Wilkinson Microwave Anisotropy Probe (WMAP) [3]. As first observed in 
1992 by the COBE satellite [8] and afterwards by several other ground-and balloon-based 
experiments, the nearly perfect black body spectrum of the CMB has little temperature 
fluctuations of the order ST/T ~ 2 72^K ~ 10~^. The angular size of these fluctuations 
encodes the density and velocity fluctuations at the surface of last scattering, with redshift 
z ~ 1100. This corresponds to the cosmological epoch when the presently observed CMB 
photons first decoupled from matter. By plotting the squared of amplitude of CMB temper- 
ature fluctuations against their wavelengths (or multipoles in an equivalent Fourier power 
spectrum), there can be allocated several peaks at different angular sizes. The position 
of the first peak is often viewed as an indicator for the spatial curvature of the universe, 
which reveals that the present universe is nearly flat and homogeneous on large cosmo- 
logical scales (> 100 Mpc), meaning that ^Itot ~ 1 with high accuracy. However, when 
assuming a flat universe only containing pressureless dust (including DM) and assuming 
the current Hubble parameter to be /i = 0.72 it 0.08 with Hq = lOOhkmsec^^ Mpc~^ (in 
agreement with observations of the Hubble Space Telescope Key project [9]), it is figured 
out that to = 9 lb IGyrs. This result, simply following from Einstein's general relativity, 
implies that a flat universe without the cosmological constant term may suffer from a seri- 
ous age problem. Introducing DE in the form of a constant A, with Oa,o — 0.73, somehow 
resolves the problem, giving to — 13.8 Gyrs with h = 0.72. 

When accepting the existence of DE, naturally the question arises, what it really is. 
Since the late 1960's when it was realized that [10] the zero point vacuum fluctuations 
in quantum field theories are Lorentz invariant, it has been attempted to associate this 
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(quantum) vacuum energy with the present value of A but without much success. Even 
when placing a cutoff at some reasonable energy scale, this quantum vacuum energy is 
still several orders of magnitude larger than the mysterious dark energy today, /9a ~ 5 x 
10~^^ GeV^ or pA ^ 10~^^^ in Planck units (for reviews, see, e.g. [11, 12]). Apparently, 
p]l^ is fifteen orders of magnitude smaller than the clectrowcak scale, m^-^r ~ 10"*^^ eV. 
No theoretical model, not even the most sophisticated, such as supersymmetry or string 
theory, is able to explain the presence of a small positive A. 

Another hurdle in understanding the nature of dark energy is that only a very small 
window in the magnitude of the cosmological constant allows the universe to develop as 
it obviously has. It is still a mystery why has the value it has today. It could have 
been several magnitudes of order larger or smaller than the matter density today, instead 
of Cl\ 2± SJljj^. This is known as cosmological coincidence problem. 

At present the most common view is that dark energy is presumably constant and 
has a constant equation of state, w^^ = — 1. But there remains the possibility that the 
cosmological constant (or the gravitational vacuum energy) is fundamentally variable. In 
a more realistic picture, at least, from field theoretic viewpoints, dark energy should be 
dynamical in nature [13]. This is the case, for instance, with all time dependent solutions 
arising out of evolving scalar fields, with an accelerated expansion coming from modified 
gravity models, holographic dark energy, and the likes. 

Interestingly enough, the recent observations (WMAP+SDSS [3]) only demand that 
— 1.04 < Wj^^ < — 0.82. In view of this wide range for the present value of dark energy 
equation of state (EoS), it is certainly worth constructing an explicit model cosmology, 
where dark energy arises because of a dynamically evolving scalar field, and see what other 
consequences would arise from such a modification of Einstein's general relativity. 

2. Constructing Inflationary Cosmology 

A complete model of the universe should perhaps feature a period of inflation in a distant 
past, leading to a generation of density (or scalar) perturbations via quantum fluctuations. 
This expectation has now received considerable observational support from measurements 
of anisotropies in the CMB as detected by WMAP and other experiments. 

In the simplest class of inflationary models, inflation is described by a single scalar (or 
an inflaton) field (?!), with some potential V{(f)). The corresponding action is 

S = l d'x^g (A _ 1 (a</,)2 _ v{<i>)^ , (2.1) 

where k = m'^ = (Stt Gn^'^ is the inverse Planck mass, with being Newton's con- 
stant, y/—g = deig^y is the determinant of the metric tensor. 

Constructing concrete models of infiation and matching them to the CMB and large 
scale structure (LSS) experiments has become one of the major pursuits in cosmology. 
Most earlier studies regarding the form of an inflationary potential relied on a prior choice 
of the potential V^cf)), or on slow-roll approximations in the calculation of power spectra 
and their relation to the mass of the field </> during infiation (see [14] for a review). The 
latter approach can at best produce the tail of an inflationary potential, but not its full 
shape [15]. Indeed, recent studies show that the type or variety of scalar potential allowed 
by array of WMAP inflationary data is still large [16]. Although, in order to understand 
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the dynamics of inflation, the idea of utiUzing one or the other form of the scalar field 
potential (motivated by physics beyond the standard model or even by theories of higher 
dimensional gravity, such as, string theory) is not bad at all, there might exist a more 
elegant way of confronting the WMAP inflationary data with a theoretical model. 

In this paper we present a different and robust approach to tackle this problem: we 
do not make a specific choice for V{4)), rather we make a simple ansatz for the scalar 
field (f) and then construct an inflationary potential, using the symmetry of Einstein's fleld 
equations. Our approach would be novel in the sense that it provides a unique shape (and 
slope) to the scalar (or inflaton) potential. The model also makes falsifiable predictions. 
The basic ideas and some of the results were presented in a recent paper [17]. 

For simplicity, we consider a spatially flat Priedmann-Robertson- Walker spacetime. 
The evolution of the fleld is then described by the equation (see, e.g. [18]) 

= mp{-2Hf/^ = -2ml-^H{4>) (2.2) 

and the evolution of the scalar potential V{(f)) is governed by 



nip acp ^ 



2 



(2.3) 



where H{(j){t)) = a/a is the Hubble parameter and a{t) is the FRW scale factor, and the 
dot denotes a derivative with respect to the cosmic time t. 

Let us first briefly discuss how the model that we are going to construct could satisfy 
inflationary constraints from the WMAP and other experiments. First, note that the term 



2ml 



is usually non- negligible (as compared to 3if^(0)) at the onset of inflation. This would 

be the case, for instance, if the mass of the inflaton field, m^ = [(fV {(/))/ dcp'^ Y^'^ , is large 

enough initially, m^ ~ mp. Once the fleld (f) rolls satisfying \(p\ <C 3H\(j)\, or equivalently, 
oc ln[a(t)], the scalar potential is well approximated by an exponential term: 

F(</)) ~ 3m^i?2(</,) ~ ^(6 - (j)'^) e2«'^' e"''^''^ ~ ^(6 - a'^)e'^''e'""^, (2.4) 

where K(j)' = a is the slope of the potential, during a slow-roll regime. The condition 
< y/6H holds in general, so V{<^) > 0. Inflation occurs as long as the condition 



a " ""P \ d(j) J 



^=H\cf>)-2ml{'^] >0 



holds, meaning that V{cl)) > (fP'. But, after a sufficient number of e-folds of expansion, infla- 
tion has to end. This is possible when the quantity {mp/ H{(p))(dH{(p) / dcp) becomes compa- 
rable to (or even larger than) unity. Recent results from WMAP [3] indicate that the spec- 
tral index of the scalar perturbations is consistent with almost flat one, = CDSSlg Q^g. 
To a good approximation, 1 — ~ a^, implying that a < 0.25. This simple picture has 
obvious and intuitive appeal, which can be realized through an explicit construction. 
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To illustrate the construction, we make the following ansatz 



= const -alnf^Vf-l'- (2-5) 



uId \ a,- / \ a- 



where a- is the initial value of the scale factor before inflation, and a and /? are free 
parameters for now. We take /3 < 0, so that after few number of e-folds, since a S> a-, the 
inflaton cf) naturally satisfies = 777^3^^ — —ot- One may think that the above choice for 
(f) is ad hoc and/or no more motivated than a particular choice of V {({)), but it is not exactly! 
Indeed (2.5) is the property of an inflaton field in many well motivated inflationary models 
that satisfy slow roll conditions, after a few e-folds of inflation. It can also be compared 
to a generic solution for a dilaton (or modulus field), i.e. (p{t) ~ + Q^Int + a^jV 
(where 7 > 0), in four-dimensional superstring models (sec, e.g. [18,19]). Additionally, the 
ansatz (2.5) allows us to construct an explicit inflationary model, providing an appropriate 
shape (and slope) to the scalar field potential. 

The evolution of as given in eq. (2.5) is provided by the Hubble parameter 



■— 7V(0) - «e-^^('^) - ^ e-2W) 



(2.6) 



where A''(^) = ln(a/aj, a = a(^(i)) and is an integration constant. We can easily 
evaluate the following two inflationary variables 



^M-^<1^,{^)-^«-JT:^ (2.8) 

(which are first-order in slow roll approximations) . The magnitude of these quantities must 
be much smaller than unity, during inflation, in order to get a sufficient number of e-folds of 
expansion, like Me = ln(a^/aj) > 50. More precisely, we require |e^| <^ 1, |r/^| < 1, except 
near to the exit from inflation where ^ 1- One may actually demand that < < 3, 
so that the scalar field potential 

V{(t>)=mlH'{cl>){3-ejj) (2.9) 

is non-negative. A typical shape of this potential is depicted in Fig. 1. The magnitude of 
(cf eq. (2.6)) can be fixed using the amplitude of density perturbations observed at the 
COBE experiments, using the normalization [20]: 

(dF/#)-^l^^/V(\/757rm|i) ~ 1.92 x 10'^. 

Typically, with a ^ 0.2 and Me = ln{ajr/a-) ^ 55, we find (assuming that (3 < 0) 

i^o ~ 7.42 X 10-^ mp. (2.10) 

This is a perfectly reasonable value, which also characterizes the average energy scale of 
inflation in most inflationary models. 
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Figure 1: The shape of the potential, for some representative values of a = 0.2, 0.3, 0.4 (top to 
bottom), P — —0.2 and N{(j)) = Ina + c. We have taken c— —10. 



As long as the parameter is slowly varying, the scalar curvature perturbation 

can be shown to be [21] 

Pi/\k) = 2^-3/24^(1 - e„)^-^/2 ( _Hl\ (2.11) 
where v = 3/2 + l/(p — 1) and a oc . The scalar spectral index for Pj^ is defined by 

The fiuctuation power spectrum is in general a function of wave number A;, and is evaluated 
when a given comoving mode crosses outside the (cosmological) horizon during inflation: 
k = aH = a^H{(l))e~^^ is, by definition, a scale matching condition and is the value 
of the scale factor at the end of inflation. Instead of specifying the fluctuation amplitude 
directly as a function of A;, it is convenient to specify it as a function of the number of 
e-folds Me of expansion between the epoch when the horizon scale modes left the horizon 
and the end of inflation. To leading order in slow roll parameters, Ug is given by [14] 

n - 1 = 2,^ - = _^!!±3^^!^i±3V±2^^^ (2.13) 

a + ^ 

where /x = Pe^^^. In the conventional case that /3 = 0, which corresponds to a scenario 
where inflation is driven by a simple exponential potential, V{(j)) obtain a 

well known result that 1 — ~ a^. Here we shall assume that Me > 47 and f3 < 0. 
Let us also define the slope or running of the spectral index n^, which is given by 

~ _ dn^ _ dus dN del) 
~ dink ~ dN dc/} dink ^ ' ^ 

(tilde is introduced here to avoid confusion with the exponent parameter a introduced in 
eq. (2.5)), where (p and k are related by 



dink ^(l-e^(0))' 
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while N and (j) are related by 

dN 

These relations hold independent of our ansatz (2.5). 

r 

.325 - 

0.275- 

0.25 
.225 



rup— = 7^^^- (2.16) 



ric 



0.955 0.965 0.97 vS^75 0.98 

. 175 ; 

Figure 2: The tensor-to-scalar ratio r ~ 16e^ vs the scalar spectral index with a — 0.2,0.18 
and 0.16 (top to bottom) and /3 = (-0.5,0). The solid (dotted) hnes are for Af^ = 60 {J\fe = 47). 

The WMAP bound on the tensor-to-scalar ratio, r ~ 16e^ < 0.3 (95% confidence 
level), implies e^j < 0.0187. This bound is satisfied for 

a < 0.1936 and > 0.9624. (2.17) 

The spectral index obtained in this way is within the range indicated by three year WMAP 
results [3] 

= 0.958l°:°J^. (2.18) 

Of course, one may directly use the above bound for and find the corresponding bound 
on r. Again by demanding that Me ^ 50 and /3 < —0.2, we find 

a = 0.2213;S, r = 0.3918;°:!^°°. (2.19) 

The smaller is the value of a, the smaller will be the tensor-to-scalar ratio (see Fig. 2), 
allowing only a small running of spectral index. For instance, if \a\ < 0.1 and /? < —0.2, 
then we find 

n, > 0.98, r < 0.08. (2.20) 

The WMAP data requires a spectral index that is significantly less than the Harrison- 
Zel'dovich-Peebles scale-invariant spectrum (n^ = 1, r = 0). Thus, given that /? < 0, 
consistency of our model (with WMAP result) seems to require 0.16 < a < 0.26. 
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Figure 3: Contour plots for — 0.95 with Afe — 50 (left plot) and Afe — 60 (right plot). 

It is also significant to note that, for a <C |/?|, there exists a small window in the 
parameter space where 



n ~ 

s 



0.95, 



r ~ 



0(10- 



10" 



in which case, however, the slope parameters a and /3 must be finely tuned. In Fig. 3 we 
show the contour plots with Me = 50 and Me = 60, representing such a case. In fact, in 
the case \a\ < 0.05, the gravity waves (or tensor modes) are almost nonexistent. On the 
right plot in Fig. 4 we show the running of spectral index 5, which is always very small in 
the parameter range 0.1 < a < 1 and /? < 0. 
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Figure 4: The scalar spectral index (left plot) and its running a with respect to /? and a = 
0.16,0.20 and 0.24 (top to bottom). The solid (dotted) lines are for Me = 60 (Me = 47). The 
running of Ug could be large only if |a| < |/3|; for example, a ~ —0.004 for a ~ 0.01 and /? ~ —0.05. 

In a model with more than one scalar field, the dependency of inflationary variables 
like and r on the slope parameters a and /? could be more complicated than the sim- 
plest explanation provided above. Nonetheless, our approach has great significance as it 
generically leads to a spectrum of primordial scalar fluctuations that is slightly red-tilted 
(n < 1) and hence compatible with WMAP inflationary data. 
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3. Constructing Quintessence Cosmology 



It is reasonable to assume that a late time acceleration of the universe is driven by the same 

mechanism usually exploited to give early universe inflation, where the potential energy 
of a scalar field dominates its kinetic term. To this end, let us assume that the current 
expansion of the universe can be described by the action 

S = 5grav + Sm = j d^Xy/^ " ^ (^^)' " ^^'5)) + 5^, (3-1) 

where Q is a fundamental scalar (or dark energy) field, V{Q) is its potential, Sgrav is the 
gravitational part of the action, Sm is the matter action describing the dynamics of ordinary 
fields (matter and radiation) and V represents a four-dimensional covariant derivative. The 
matter part of the action (3.1) can be written as 

Sm = j d^xC{'il:^,A\Q)g^,) = j d^x^A\Q) ^p^, (3.2) 

where ipm represents collectively the matter degrees of freedom and radiation. In the 
above definition of the matter Lagrangian, the implicit assumption is that matter couples 
to g^i, = A{Q)'^g^v, rather than the Einstein metric g^n alone. This assumption then results 
in a non-minimal coupling between the scalar field Q and matter components {pi). The 
matter-scalar coupling A{Q) may be understood as a natural modification of Einstein's 
GR which can be motivated by, for instance, scalar-tensor theory. For further discussions 
on theoretical motivations of this coupling, see, for example, [23-26] . 
The coupling A{Q) actually generates a new term, namely 

2 5C!i _ dA{Q) 



^ 5Q dQ 

in the scalar wave equation for Q. This expression also implies that radiation does not 
couple to the scalar field Q since its trace of the energy-momentum tensor equals zero. As 
we will show the coupling ^^^^ introduces several qualitatively new cosmological features. 

As is well known, the cosmological constant case (or more generally Einstein gravity 
with a cosmological term) arises as a special limit of the present model, for which 

and hence V{Q) = const = A and A{Q) = const. The model then reduces to the ACDM 
cosmology, given that dark matter is characterized by non-relativistic particles alone, = 
0. The cosmological term A, which is governed by the equation 

mlG'^ = 2A6i^ + Tlf, (3.5) 

can clearly act as a source of gravitational repulsion or putative dark energy. 

All the discussions so far have been made without making any particular choice of 
metric. Thus the nature of V{Q) acting as a repulsive force is rather general. For a more 
detailed treatment, it is necessary to evaluate the equations generated by variation of the 
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total action S = S„^^ + . Therefore a particular choice of a metric has to be made. We 

grdyV Til 

make rather standard choice of a spatially flat FRW metric: 

ds'^ = -dt^ + a (tf dx^, (3.6) 

where a{t) is the scale factor of a FRW universe. This choice of the line element is well 
motivated by the observational fact that the universe is spatially flat on largest scales, which 
is consistent with the concept of inflation, discussed in the previous section. Of course, this 
choice of metric may lead to systematic errors in the calculation, as the universe actually is 
not homogenous at smaller (or galactic) scales, as pointed out, for example, in [22], which 
is ignored in this simplified assumption. 

In the minimal coupling case, A{Q) = 1, it is easy to see that 

p,cx[a(t)]-3(i+-'), (3.7) 

where w- = Pi/pi- In the non-minimal coupling case the modified scale factor a is given 
by a = a{t)A{Q). As a consequence, different equation of state parameters (cf eq. (3.3)) 
would cause different energy densities to evolve differently with changing scale factor: 

p,cx{a{t)A{Q))-^('+^i). (3.8) 

This implies that pm oc {a{t)A{Q))~'^ and pr oc {a{t)A{Q))~^, respectively, for ordinary 
matter and radiation. It also shows that radiation never directly couples to the scalar field, 
even with A(Q) being an arbitrary function of Q. As explained in [17], the coupling A(Q) 
can be relevant, especially, in a background where p^ is much larger than p^^^^ (where 
'^crit = 3i?^/87rG), e.g., a galactic environment. 

3.1 Basic Equations 

Taking a variation of the action (3.1) with respect to g^^'^ and then evaluating the tt and 
XX components of Einstein's equation leads to the following two equations (cf eq. (A.l)): 

-^H^ + Iq^ + V (Q) + A^ (Q) ^Pi = (3.9) 

% 

^2 (2 + 3 + ^ g2 _ y (g) + ^4 (g) ^ ^ (3 -^0) 

i 

A variation with respect to the scalar field Q, while considering an explicit matter-scalar 
coupling, yields the following equation of motion for Q (cf eq. (A. 2)): 

i 

the so-called the Klein-Gordon equation for Q. It shows that the scalar field Q couples to 
the trace of the energy-momentum tensor g^^T^^ satisfying 

-V^Q = Q + ^HQ = - A^ (3-12) 

dQ dQ 1^ 

There is dissension about the sign of the coupling term between the scalar field and matter 
in above equation in the way that it might be +A? '^^^q^ instead of —A? '^^q^ T^. In 
this paper, the negative sign, as written in eq. (3.11), will be used. 
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The above set of equations can be supplemented by a fourth equation, arising from 
the equation of motion for a perfect barotropic fluid 

This finally leads to (cf eq. (A.7), see also [27]) 

p^ + 3i/(l + u;Jp^ = ^^(l-3u,Jp, (3.14) 

Out of the four equations (3.9)-(3.11) and (3.14), only three are independent, meaning the 
conservation equation of the perfect fluid (3.14) can be derived without the assumption 
(3.13) but only by combining (3.9)-(3.11). General covariance requires the conservation of 
the total energy density, p^^^ = Pq + ^Pi, which is obviously the case in our model (see 
also the appendix in [28]). 

Next we make the following substitutions: 

(3.15, 

and 

_ \Q^-v{Q) _PQ ..... 

Wr, = = — • (3.16) 

\Q'' + V{Q) PQ 

These substitutions and further simplifications lead to the set of four equations: 

Y,^^ + ^Q = l (3.17) 

i 

2e + 3 (l + u;(j) Oq + 3^(1 + u;.) = (3.18) 

i 

Q.'Q + 2eVtQ + mQ (i + w;q) ^(r?iQ.) =0 (3.19) 

i 

^ f); + 2 6 ^ f], + 3 ^ (J^, (1 + t«J ) - Q' 5^ (r/, J7 .) = 0, (3.20) 

i i i i 

where, as above, the prime denotes a derivative with respect to e-folding time N = ln[a(t)] + 
const, Q' = Q/H, rj- = 1 — Sw^^ and ag = In the above we have used the relation 

did , , 

Equations (3.17)-(3.20) represent the most general case of an evolving universe based on 

the general action (3.1). Equation (3.17) is simply the Friedmann constraint for the as- 
sumed flat universe. Changing the sign of the coupling Oq to the trace of the energy- 
momentum tensor in eq. (3.12) would cause a change of sign from -j-Q' ag (r/i Q,-^ 
to -Q'ag Y.i iVi^i) in eq. (3.19). Adding eqs. (3.19) and (3.20), we find 
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which can be interpreted as a global energy conservation equation. Thus, for not violating 
this principle of energy conservation (3.22), a sign change in (3.19) automatically implies 
a change in (3.20) as well. 

When having a particular solution of the equations (3.17)-(3.20), it is of great interest 
to study how the corresponding potential looks like and how it affects the cosmic evolution 
of our universe. From the last expression in eq. (3.15), we find 

- tt2 f^^Q 1 ^r,'\^ 



K 



V{Q)^H^ l^- {Q'Y), (3.23) 



which will be used later. Of course, in the case of a minimal coupling {A{Q) = 1), ag 
vanishes, reducing the number of degrees of freedom in the system of equations (3.17)-(3.20) 
by one, which then makes the system easier to handle. Anyhow in both cases (ag = and 
Uq 7^ 0) it is not possible to find an analytical solution of this system without making some 
additional assumptions as there are more degrees of freedom than independent equations. 
In fact, the number of degrees of freedom depends on the number of matter components 
included in the analysis. 

As the first check for compatibility of the model, it is useful to consider some simplified 
solution of the equations (3. 14)- (3. 20), by expressing all matter fields as one component, 
w- = w . By applying eq. (3.21) to eq. (3.14), and after a simple integration, we get 



(0) g-31na 



J (<3'«Q {^w^ - 1) dlna 



(3.24) 



with p'^ being an arbitrary constant. The coupling Og may be constrained by observations 
perhaps only in the combination Q'cxq. One can study the effect of this coupling on both 
CMB temperature anisotropics and evolution of linear matter perturbations, as in [29]. In 
the minimal coupling case, one has 



-3(l+«; 



(3.25) 



This is exactly the behaviour one would expect from general relativity. Equation (3.25) 
yields oc a"^ in a universe containing only ordinary matter (or dust), while for radiation 
Pr OC a~^. Transposing cq. (3.18) leads to a general expression for the equation of state of 
the DE component which can generally be written as 

where all possible forms of matter are included, e.g. pressureless dust {w-ai = 0), radiation 

(wr = 1/3), stiff matter (tfsm = 1), domain walls (wdw = ^2/3), etc. 

For further analysis, it is useful to introduce the so called effective equation of state 
parameter w^^, which is defined by 

w^eff = — ' Ptot =PQ+Pi, Ptot = PQ + Pi, (3.27) 

Ptot 

whereas pQ and pQ are as defined in eq. (3.16), and p- = A'^{Q)p., p. = A'^{Q)p,-. The mean- 
ing of u)gg is somehow that of a mean equation of state of all matter-energy components. 
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including the dark energy component pg. Prom eq. (3.18), together with eqs.(3.15)-(3.16), 
we find that the total pressure is given by 



ftot = ^ i-^-^Q-E^i]- (3-28) 



Combining this expression for p^.^^ with the expression of total energy density p^^^, as defined 
by (3.27), and using again the substitutions (3.15)-(3.16), we get 



^eff = -1 - y • (3.29) 



This expression is valid even in the most general non-minimal coupling case. Similarly, we 
can define the deceleration parameter q as 

q = ^ = -!-€, (3.30) 

showing that the name "deceleration parameter" makes sense in such a way that q > for 

a decelerated expansion (a < 0), while g' < for an accelerating expansion (a > 0). 

As explained in the Introduction, one reason for considering a universe containing a 
nonzero DE component, either in the form of a cosmological constant A or a dynamically 
evolving scalar field Q, is the recently observed accelerated expansion of the universe by the 
Supernova Cosmology Project and the High-redshift Supernova Search team [1]. Thus we 
find it useful to study the solution of equations (3.17) to (3.20), which yields an accelerated 
expansion in a general context of nontrivial matter-scalar coupling. 

Indeed, independent of any assumption or specific composition of the universe, simply 
the condition w^g < —1/3 at some stage of cosmic evolution yields an accelerated solution. 
It should be noted at this stage, that no such general connection can be established between 
the DE equation of state parameter Wq having a specific value (even like Wq = — 1) and 
the universe being in an accelerating phase. 

In the notations used in this paper, both the non-baryonic (cold) dark matter and 
ordinary matter (prcssurclcss dust) are combined in one matter constituent Jl^^^. As w^,^^ ~ 
is a rather good approximation for the equation of state of cold dark matter (since it is 
non-relativistic) this combination seems to be reasonable. This assumption as regards the 
composition of the universe today implies that its only constituents are cold dark matter, 
ordinary matter, radiation and DE. Putting this composition (0^^^ = 0.27, 0,q = 0.73 and 
= 10~^) of today's universe into the very general expression of the DE equation of state 
parameter Wq (cf eq. (3.26)) and using again = and Wq ~ —1, the value e = —0.4 
is obtained. This implies that for Wq at least being close to —1 the universe is in an 
accelerating phase today, which is what is observed. 

The general considerations so far seem to be consistent with observations. As observa- 
tions seem to indicate a value for Wq close to —1, the possibility of a dark energy component 
simply being a cosmological constant cannot be ruled out. But it is also important to re- 
alize that the effects of a slowly rolling scalar field would be almost indistinguishable from 
that of a pure cosmological constant if kQ' = m~^{Q/H) < 0.1 at present. Evidence 
for Wq ~ — 1 could actually imply that the field Q is rolling only with a tiny velocity at 
present. This point should be more clear from the discussion below. 

All the examinations so far have been in a rather general way without imposing any 
additional assumptions. For sure that is not really satisfying, as one might be interested in 
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an analytic solution of the system of equations (3.17) to (3.20). As mentioned above this is 
not possible without further input because of the number of degrees of freedom exceeding 
the number of independent equations. In the next two subsections two different analytic 
solutions will be presented making some simple additional assumptions. According to the 
present constitution of the universe being ~ 0.27 and 0,q ~ 0.73, it is reasonable to 
neglect the radiation component at least for redshift z < 0(10). Therefore the model 
universe assumed in the next two sections is thought to only consist of cold dark matter 
and ordinary matter combined in one component with a common equation of state Wm = 
and a DE component represented by the scalar field Q with a variable EoS Wq. 
The system of equations (3.17)-(3.20) can then be expressed in the form: 

^^„, + nQ = i (3.31) 

2e + 3 (i + Wq^ nQ + 3n^ = (3.32) 

n'Q + 2enQ + 3nQ (i + wqJ +Q'aQ^^ = o (3.33) 

^L. + 2 e + 3 - Q' = 0. (3.34) 

The number of free parameters in this system is five {'^^, ^q, Wq, e and cCq), meaning two 
additional assumptions have to be made to find an analytic solution. To proceed further, 
we make the following assumption: 

Q = Qo + mp a ln[a(i)] = Qq + rupa {N + const) , (3.35) 

where a is a constant which needs to be fixed by observations. This relation actually 
represents a generic situation that the field Q is rolling with a constant velocity, Q' = 
const ^. In the minimal coupling case this is enough, while in the non-minimal case (a^ / 0) 
one more assumption is required, which will be discussed below. 

Simply transposing (3.16) and utilizing the relation between ^ and as given by 
(3.21), yields the following useful relation 

^^g^-sn^^cl-sn^ (3.36) 

Supplementing equations (3.31)-(3.34) with this equation is an elegant way of imposing an 
additional constraint into the model. 

3.2 Uncoupled Quintessence 

In the A{Q) = 1 case, the system of equations (3.31)-(3.34), supplemented by eq. (3.35), 
can now be solved analytically. The explicit solution is given by 

= ^ " c,Aexp[AiV] + 3 ^^"^^^ 
- -c,A^ 

"'Q " 3a2exp[-AiV] + 3Aci ^^'^^^ 
_ «^ A - 9 exp [-\N] 
^ 2c,A + 6exp[-AiV] ' ^ ' 



'^Note that we are demanding Q' = dQ/d\na = Q/H ~ const, not Q = const. For a < 0.6, our approach 
appears to give consistent results when applied to observational data; see also the review [30] for extensive 
discussions on various methods of reconstructing dark energy potentials. See ref. [31] for a very different 
approach of dark energy reconstruction. 
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where A'' = N{Q) = ln[a{Q{t))] and we have made the substitution 



A = 3 - a^. (3.40) 

Using eqs. (3.29) and (3.30) we also evaluate 

_ 3 exp [-AA^] - (6 - 5 + q,4) 



(3.41) 



6 exp [-XN] + 2 A 
9 exp[-AiVf+3c^A' 



This general solution contains three free parameters (A'", a and c^). To keep the solution 
as general as possible it is useful to just fix one free parameter in terms of the other two. 
The integration constant can be fixed in terms of the field velocity a by using the 
observational input = 0.27 at present. The e-folding time N in relation to the cosmic 
time t is only defined up to an arbitrary constant, so it needs to be normalised in some 
way. For simplicity, this will be done by taking A'^ = at present. Thus, the condition 
^JN = 0,a,c^] = n^Q yields 

which now makes it possible to express eqs. (3.37)-(3.39) just in terms of the two free 

parameters A'^ and a. For further analysis it is useful to parameterize the solution in terms 
of redshift z. By utilising the dependence of the redshift on the scale factor a(t), it is easy 
to obtain the relation between A^ and z: 

^^^^A^^ ao^ exp[Ar] = ^ = ^. (3.44) 

Aem a{t) ao l + z 

Here t is the time when light was emitted, that is observed now. Thus choosing ag to be 
the present scale factor automatically implies the normalization A'^ = at z = 0. From 
eqs. (3.37)-(3.39), one can clearly see that the solution is symmetric in a as only even 
powers of a occur. Thus, without loss of generality, in the further analysis only positive a 
will be considered. For the value of c-^ satisfying (3.43), and fi^Q ~ 0.27, it is easy to see 
that = for = 2.19, implying that a^^-^^ = 1.48. This value of a has some significance, 
when looking at the evolution of Wq for different values of a. 

From the left plot in Fig. 5 it is easily seen that a = 1.48 yields Wq = whereas Wq < 
for a < 1.48. Here a = 1.48 separates solutions with the energy described by means of 
Q being attractive or repulsive. Thus, if il^ = 0.27 at z = 0, then a solution with self- 
repulsive DE requires < a < 1.48, whereas a = equals the cosmological constant case 
with Wq = —1, which can be also seen in Fig. 5. WMAP data combined with the Supernova 
Legacy Survey (SNLS) data yields a significant constraint on the equation of state of the 
dark energy, —1.04 < Wq < —0.82 (with 95% CL), sec also refs. [32,33]. However, here 

we consider only the region — 1 < Wq < —0.82 so that > 0, that is, without going 
to a phantom regime. This would require a to be in the interval < a < 0.63 which 
can also be inferred from Fig. 5. What can also be seen in Fig. 5 is that for all a in 
the range < a < 1.48, Wq for z oo, thus implying the DE component being 
indistinguishable from pressureless dust for high redshifts and only becoming the observed 
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Figure 7: The effective EoS w^^ with respect to redshift z: (left plot) a — 0, 0.7, 1.0, 1.48, 1.6 
(bottom to top) and (right plot) a = Q, 0.4, 0.6, 0.8, 1.0 (bottom to top). 



self-repulsive form of energy in recent time. For a further understanding of this solution it 
is useful to look at the deceleration parameter q. 

In Fig. 6, it can be easily seen that, for all < a < 1.48, q gets negative somewhen 
between redshifts z = and z = 1, which implies that in this model accelerated expansion 
is a rather late time phenomenon with the universe getting into an accelerated phase the 
earliest for a = 0, corresponding to the cosmological constant case. In the case of a = 1.48, 
q exactly equals 0.5, corresponding to a decelerated expansion at constant deceleration. 
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Figure 8: il^ and fig with respect to a and z. These quantities may not change with z only if 
a = Q;^j.j^ = 1.48, in which case obviously there won't be a cosmic acceleration. 

Finally, for a > a^^,-^ = 1.48, q is greater than 0.5 and increases with decreasing redshift, 
yielding a decelerated expansion. This fits to the evolution of the dark energy EoS Wq, as 
seen in Fig. 5 (for a > a^j.^, Wq > 0). 

From Figs. 7 and 8 we can see that for the solution which leads to a late time ac- 
celeration {w^g < —1/3) the universe is clearly dominated at high redshift by with a 
transition to Qq dominance in recent time leading to = 0.27 and 0,q = 0.73 at z = 0. 
(That for sure does not come surprisingly, since that was the assumption made when fixing 
c^). It is perhaps more interesting to note that for a = 1.48 the ratio fig/ri^ remains 
constant for all z, whereas, for a < 1.48, the early universe would be dominated by 
with a shift to dark energy dominance in the recent epoch. The observed acceleration and 
DE dominance correspond best to values of a closer to zero. 

Uncertainties in the current value of 0^ affect a^it) to some extent, and hence the 
predicted value of Wq at some fixed redshift. That is, for a value of $7^ different from 0.27 

at present, the critical value of a, i.e. a can also be different. However, the 

general bahaviour of the solution would be similar. 




Figure 9: ilq with respect to and q, for a varying z ~ {5,0} and a = {1,0} (left plot) and 
a = 0, 0.15, 0.3, 0.45, 0.6 (right plot, bottom to top), lowers to —1 at a low redshift. 

The left plot in Fig. 9 is a three-dimensional illustration of the above discussed fact. 
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that a transition to the accelerated phase {q < 0) occurs for Wq tending to —1 and fig 
tending to +1. The right plot in Fig. 9 is a two-dimensional projection of the latter and 
thus just gives another illustration of the already discussed relation between Wq and Qq 
for the accelerating case, where only accelerating solutions with a < 0.6, which actually 
lead to Wq < —0.83 at z = 0, are examined. 



V 
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Figure 10: (Left plot) 



'vjQ) 



with respect to z and a. (Right plot) 



'vjQ) 



0, 0.4, 0.8 and 1.2 (from bottom to top, at the right end of the graph). 



with respect to z for 



The discussion so far has been based on the idea of a dark energy as described by 
the scalar field Q with some potential V{Q). For obtaining the analytical solution, (3.37)- 
(3.39), no particular choice was made for the potential. The only one assumption made 
was that the field might be rolling with a constant velocity a, with respect to the e-folding 
time = In a. Thus it would be worth looking at the shape of the potential as determined 
by this particular solution, following the idea of reconstruction underlying the focus of this 
paper. For obtaining the analytic expression of V{Q), it is useful to consider the set of 
substitutions made in (3.15)-(3.16). By utilising the additional constraint (3.36), it is easy 
to see that 



Y = 



V{Q) 



3 a 



a 

T 



Y is actually a dimensionless variable, which takes the value y = 3 in a pure de Sitter 
space. The variation of Y shown in Fig. 10 seems quite natural and can be understood in 
the following way. In order to get an accelerated expansion of the universe, with Wq close 
to —1 at a low redshift, Y/3 should exceed fi^ in the recent past. 

In order to find the potential, it is necessary first to evaluate the Hubble parameter 
H, which can be easily done by solving the equation 

e[N]H[N]= H'[N]. (3.46) 

The analytic expression of H is given by 



H 



^2 exp 



3 exp [-\N] + A, 



(3.47) 



The numerical constant can be fixed by the assumption that H[N = Q] = Hq. Hence 



V^3 + A ■ 



(3.48) 
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Finally, the quintessence potential takes the form 

V[Q{N)) = \cl exp [- a^N] (Sa^ exp [-\N\ + (l8 - 9 + a^)) 



(3.49) 




In Fig. 11 it is clearly seen that V{Q) increases exponentially with increasing redshift 
z, whereas this increment is more steep for larger values of a. As it should be, in the 
a = case, the potential takes a constant value. In fact, the assumption of Q rolling 
with a constant velocity {Q' = a) yields that the potential V{Q{z)) must take a shape to 
cause this behaviour for Q. An exponential shape for the potential is no surprise. The 
quintessence potential constructed in this way takes the following form 



K^V 



- o(c3+kQ) 



2>a exp 



a 



3) (c3 + kQ) 



a 



+ (18 - 9a^ + 



2 ' a^l 



(3.50) 



The integration constant Cg can be set to zero, without loss of generality, while and C2 
can be fixed in terms of a (and Hq), using eqs. (3.43) and (3.48). The potential can be 
brought into a form where it only depends on Q and a: 



V{Q) 




K exp 



[a 



3)g 



anil 



(3.51) 



This potential is clearly double exponential in form and it would find interesting appli- 
cations even for the early universe. As discussed in [34], one may be required to have 
a < 0.8 in order to satisfy the bound on during big bang nucleosynthesis, namely 
Oq(1 MeV) < 0.1. It is also interesting to note that such a potential can easily arise from 
some fundamental theories of gravity in higher dimensions (see e.g. [35]). 

The form of the potential, as it can be seen in Fig. 12, is not surprising, as it allows 
the field to "roll down" the slope of a decreasing V{Q) for an increasing Q, or a decreasing 
redshift. As expected, the slope of the potential is shallower for smaller values of a and 
equals zero in the cosmological constant case, a = 0. 
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with 



respect to kQ and Cg, for a ~ 0.6. 



3.3 Coupled Quintessence 

As already mentioned above, when solving the system of equations (3.31)-(3.34) with the 
additional constraint (3.36) in the general case (ag 7^ 0), one more constraint is needed 
to get an analytic solution. It is most canonical to assume k^cHq = const = X) which 
represents the case of so-called exponential coupling between the scalar field Q and matter, 
as A(Q) tx e^^'^/''^p\ This additional assumption then leads to a general analytic solution 



C 



3 + C, (3 + QX 



q2) e' 



w, 



c, C exp [QN] ia^ - 3) + 3 xa 



Q 



3C4 C exp [CN] + 3a2 
-3C 



3x" 



2 2C4C exp[CAf] +6' 

where 

C = 3 + xa — 

Further, the analytic expressions for q and w^q are given by 

_ . «\ 3C 
^ 2 ^ 2C4C exp [CiV] + 6 



w 



eff 



1 ^ 2C 

' ^ 3 ^ 2C4C exp[CiV] +6' 



By solving the differential equation (3.46), the Hubble parameter is found to be 

■A^(3 + ax)" 



H{Q) = Cg exp 



C4 C exp [QN] + 3, 



(3.52) 
(3.53) 
(3.54) 

(3.55) 

(3.56) 
(3.57) 

(3.58) 



where N = N{Q). The integration constant can be fixed by the assumption that 
H[N = G\=Hq. This yields 
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One normalizes N such that = corresponds to a = Bq = 1. Further, insisting that 
(^N = 0, a, c^,x) = at z = fixes the integration constant in terms of a and x- 

Compared to the minimal coupling case (x = 0), now the symmetry in the solution between 
positive and negative a is lost. However, a simultaneous change in sign of the parameters 
a and x keeps unchanged. Thus in further analysis only the properties of a solution 
with positive a but either sign of x will be examined. In the discussion that follows the 
case X < will characterize solutions with a and x having the opposite sign, while the case 
X > will characterize solutions with a and x having the same sign. 

As the parameters q, w^f^ and e are all intimately connected by eqs. (3.29) and (3.30), 
only the x-dependence of dark energy EoS Wq will be examined as an exemplary. 

As can be seen from Fig. 13, the x > solution decreases Wq, whereas the x < 
solution increases Wq (for fixed a and z). That is, a negative x causes the universe to get 
into an accelerating phase later than for positive x- In general, the decrease in Wq would 
be steeper for x < than for x ^ 0. It is also important to realize that the coupling x 
does not affect the value of Wq at z ~ but only at higher redshifts. 




Figure 13: The dark energy equation of state Wq with respect to z and a for x = —0.4 (left plot) 
and X = +0.4 (right plot). The x > solution yields a more negative Wq at a given rcdshift. 




Figure 14: ilg with respect to z and a for x ~ —0.6 (left plot) and x = +0.6 (right plot). 

In Fig. 14 we show the variation of dark energy density with the field velocity a and 
the redshift z. It is found that, for fixed a (< a^j^), can be smaller (larger) at higher 
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redshifts for X > (x < 0)- This behavior would be somewhat opposite in an decelerating 
universe with a > a^^^^. This behaviour is expected by the x-dependence of q, since an 
increase in matter density also increases q and vice versa. For a better understanding of 
this situation, it is useful to study the behaviour of the potential V{Q). 

It is also worth examining the values of dark energy EoS Wq with a varying x- In the 
case X < 0, an increasing negative x decreases Wq, whereas an increasing positive x will 
increase Wq with respect to the value it has in the minimal coupling case, x = 0; one may 
compare the figure 15 with 5. 
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Figure 15: The dark energy EoS with respect to z. Left plot: x = +0.1 and a = 0, 0.2, 0.4, 
1.0, 1.4 (bottom to top). Right plot: a = 0.4 and x = -0-2, -0.1, 0, 0.1, 0.2 (top to bottom). 

In analogy to the previous section '^^^q^ can be obtained by using eq. (3.45). In the 
X 7^ case, the effective potential consists of V{Q) and an additional term depending on 
the matter-quintessence coupling Uq. The functional form of V^^{Q) can be obtained by 
integrating the right hand side of eq. (3.12) with respect to Q. The result is given by 



K^VJQ) = Q.bcl exp 



(3 + ax) (cg + KQ) 



a 



X C4 C exp 



C (Cq + K Q) 



a 



(e-a^) -3 (6 + 4ax-3 




(3.61) 



where Cg is an integration constant. One can fix and using eqs. (3.60) and (3.59), and 
also eq. (3.44). We exhibit the shape of this potential in Fig. 16. 

As can be easily seen in Fig. 16, for a small a (< aj.j,jf), V^q{Q) increases with an 
increasing z, which allows the field to "roll down" with a constant velocity a, with the 
slope being zero for a = 0, as in the x = case. For large a (like a > \/2), instead, V^^^iQ) 
decreases with an increasing z. This should not come as a surprise; this behaviour has 
its origin in the value of oi^^-^^ which is lowered for x < 0. For a given a, the slope of the 
potential is shallower for x > than for x < 0, with vanishing difference at lower redshifts. 

We conclude this section with the following two remarks. Firstly, in our model, it is 
possible that the current acceleration of the universe is only transient. This can easily hap- 
pen, for Qg < 0, when J Oq(1 — Stt^j)/)- dQ (where oc a~'^~^"'«) becomes comparable 
to (or exceeds) k'^V{Q), making the effective potential almost vanishing (or negative). 

Secondly, in the case both the ordinary and dark matter have same coupling with the 
quintessence field Q, current observational constraints (from Cassini experiments and the 
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Figure 16: The effective potential V^g{Q) with respect to redshift z and the slope parameter a, 
in the units = 1 = k, for x = —0.5 (left plot) and x = +0.5 (right plot). We have taken Cg = 0. 

likes) only demand that < 10~^, while this bound is significantly relaxed if dark matter 
can have much stronger coupling with Q. It should be the astrophysical observations that 
decide whether < or Og > 0. The answer to this question can have interesting 
cosmological effects which we aim to study in future work. 



4. Confronting models with data 

In this section we confront our models with recent cosmological datasets (Supernova Legacy 
Survey (SNLS) and SNIa Gold06 datasets) following the methods discussed, for example, 
in refs. [36, 37]. 

In the minimal coupling case, since pg + 3H{1 + Wq)Pq = (i.e. pq and p^ are 
separately conserved), we get 



{l + w{z^)) 



10 1 + ^1 

Without any prior on w{zi) or pg, it can be shown that 

fln{l+2) 



dzi 



(4.1) 



H{z) = Hjn^^il + zy + n exp 



and 



(1 + ^) 



{I + w{z^))dln{l + z^] 

d\nH 



1/2 



dz 



1 



(4.2) 



(4.3) 



In our model we have assumed that nipQ/H = a. In this particular case, with = 0, 
the Hubble parameter H{z) as a function of the redshift z is given by (cf eq. (3.47)) 



H{z) = Hjn^^,{i + z)3 + (1 - n^jii + .) 



(4.4) 



where J7^q = 3/(3 + c-,^). Using this expression of H(z), we show in Fig. 17 the best fit 



form of w{z) for the SNLS data with a prior 
state Wq{z) is given by 



Wq{z) 



(l-^^™,o)(«'-3) 



3(1-0 J + a20 (1 + ^ 



0.24. The dark energy equation of 

(4.5) 
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Clearly, knowledge of and a would suffice to determine Wq{z). In the a = case, 
Wq{z) = = — 1. In tables 1 and 2 we present the best fit values of a and Wq for different 
choices of O^n- 



W ( Z ) 




0.2 0.4 0.6 0.8 1 

Figure 17: The best fit form of w{z) for the SNLS datasets for a prior of il^^ — 0.24 along with the 
la errors (shaded region). The (black) solid line corresponds to the ansatz w{z^) = Wq-\-w^z^/ (X+z^) 
(cf eq. (4.3)). The three other lines correspond to a = 0.4,0.2109,0 (top to bottom) and Wq{z) 
given by eq. (4.5). With n^^^^ = 0.24, a = 0.2109 minimizes the x^„^ (= 104.18). The SNLS data 
may favour a lower value of (as compared to the Gold SNIa dataset). Further, with a canonical 
quintessence, so that Wq{z = 0) > — 1, we may require Sl^^„ < 0.2592. 



Table 1: The best fit values of Wq{z) and a for the SNLS datasets for a given Vt 



^mO 


a 


Wq{z = 0) 




0.22 


0.2948 


-0.9628 


104.23 


0.23 


0.2573 


-0.9713 


104.21 


0.24 


0.2109 


-0.9805 


104.18 


0.25 


0.1476 


-0.9903 


104.16 


0.259173 


0.014 


-0.9999 


104.14 



Table 2: The best fit values of Wq{z) and a for the Gold SNIa dataset for a given Q, 





a 


Wq{z = 0) 




0.23 


0.4001 


-0.9307 


178.64 


0.25 


0.3376 


-0.9493 


178.21 


0.27 


0.2544 


-0.9704 


177.76 


0.29 


0.0933 


-0.9959 


177.31 


0.2929 


0.0100 


-0.9999 


177.25 
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Figure 18: The best fit form of w(z) for the Gold SNIa dataset for a prior of 51^^ ~ 0.27 along with 
the la errors (shaded region) with H{z) given by eq. (4.7) (left plot) and H^^^{z) given by (4.10) 
(right plot); is minimized for a — 0.4735 and a^^ = 0.0633. The (black) solid line corresponds 
to the best fit line with x^^^i— 177) and the three other lines represent Wq{z) (cf eq. (4.8)) with 
a = 0.6, 0.4, 0.2 (top to bottom) and = x = 0.4. 



The Gold SNIa datasets could actually fit better with coupled quintessence (or interacting 
dark energy) models (cf Fig. 18). 

In the non-minimal coupling case, Pq is not separately conserved, since pg + 3H{1 + 
Wq)Pq = aqHQ' p^; of course, the total energy is always conserved: /Otot~'"'^-^('^tot"'"^'tot) ~ 
0, where p^.^^. — Pm~^ Pq- Using the relations d/dt = H{d / d\na) and Ina = — ln(l + z), 
we get 



Pq = exp 



l + Zl 



dzi 



+ 



Q'a. 



l + zi 



exp 



dzi 



dzi 



(4.6) 



In particular, with nipQ' = a and Qq 
to be 



d In A(Q) 
d{KQ) 



X, the Hubble parameter H{z) is found 

(4.7) 



where „ = 3/(3 + C4C) ^^id C = 3 — + ax- The dark energy equation of state is 



Wq{z) 



«X^n^O + (l-^mo)(«'-3)(l + ^ 

?,{l-0){l + z)-^ + a{a-x)^ 



(4. 



mO 



Next we briefly discuss about an interesting possibility (leaving the details and further 
generalization to a forthcoming paper). In the non-minimal coupling case, the Hubble 
expansion parameter that one measures (in a physical Jordan frame) could actually be 
different than the one given by (4.7) by a conformal factor. Given that 



^obs(^) 

H{z) 



exp[x(Q/?7ip)] oc exp[xa; In a] = a"^ = {1 + z) 



-"X 



(4.9) 
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we find 

^obs(^) = ^oV^™o(l + + (1 - ^mo)i^ + ^)"'^"^- (4-10) 

Using this expression of H{z), we liave presented in table 3 tlie best fit values of a and 
which minimize the for the Gold SNIa, SNIa+CMB-shift (WMAP)+ SDSS data 
sets for a given „ = 0.27. 



Table 3: The best fit values of a and q;qq, with la errors for Wq = Wq{z = 0). 





a 


«0o 


Wq, (eq. (4.8)) 


Wq, (eq. (4.3)) 


SNIa 

SNIa+WMAP+SDSS 


0.4735 
0.5142 


0.0633 
0.0583 


-0.90t°;i 
-0.88i°i? 


-0.94l°:l° 

-o.92l°:°3^ 



The mean value of Wq, obtained above is within the range indicated by WMAP3+SDSS 
observations: w-^^ = — 0.94l;toioi [3]- The best fit value of Uq is found to be Uq ~ 0.06, 
but in our model it may contain significant numerical errors, namely Uq = 0.06ib0.35, which 
thereby implies the consistency of our model with general relativity (for which Uq = 0) at 
lo" level. To illustrate this result we show in Fig. 19 the best fit plot with Uq = 0. 

w(z) 

1 , , , , , , 



0.5 




-1.5 



-2 ; ^ ^ ^ ^ ^ 

0.25 0.5 0.75 1 1.25 1.5 



Figure 19: As in Fig. 18 (right plot) but with = 0. 
The post-Newtonian parameter 7 is related to (= x) through the relation [39] 




1-7 
1 + 7' 



With the best fit value ~ 0.06, this yields I7 — 1| — 7.1 x 10 ^, which is not far from a 
constraint coming from the Solar-system experiments, i.e., I7 — 1| < 2 x 10^^. Moreover, in 
the non-minimal coupling case, with A{Q) = e'^^'^/mp) ^-^ ^ ^j^^i Q{t) = alna + const, 
there arise constraints on the time variation of Newton's constant. With a scalar field Q 
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Ggff /G 




z 



Figure 20: The time variation of Newton's constant in the non-minimal case. 

conformally coupled to the matter, the effective Newton's constant (measured, e.g., in a 
Cavendish type experiment) can be given by 

^ = A{Qf (l + al) = (1 + z)-"^(l + x'). (4.11) 

The time derivative of Newton's constant generally depends on the coupling A{Q) and 
its derivative, Og. In our model, with ax > 0, the case of decreasing Wq{z) (at a lower 
redshift) corresponds to an increasing Newton's constant that boosts cosmic acceleration. 
For the SNIa best fit value (0,0^) = (0.4735,0.0633), the variation of G^^ in the 

redshift range z = {0,20} is less than 10% (cf Fig. 20) and l^^^l/G^g = 0.029/ii7o ~ 
2.1 X 10~^'^yr~^. We should mention that the current solar system constraint on G^^/G^^ 
could be more stringent than this, namely {dG / dt) / G < 10~^^ yr""^ (see, e.g. ref. [40] 
which derives constraints on G/G and G/G for a model where Q-field is explicitly coupled to 
the Einstein-Hilbert term); it is because the relevant background when studying the solar 
system is not the cosmological but the solution of (3.12) corresponding to the galactic 
environment, where Q/H ^ and p^^^ » p^^,-^ = SH^/SirG. In order to properly address 
the question of time derivative (or variation) of Newton's constant, one has to consider in 
detail the dynamical system where Oq is time-varying. This is left for future studies. 



5. Conclusion 

In this paper we have outlined construction of an effective cosmological model each for 
inflation and dark energy (or quintessence), within the framework of the standard scalar- 
tensor theory. The general assumption has been that the evolution of our universe can be 
described by Einstein's gravity coupled to a fundamental scalar field plus matter, described 
by the general action (3.1). The gravitational part of the action, which is important for 
constructing a model of inflation, contains a scalar field lagrangian. The matter part of 
the action contains all possible matter constituents in the form of a perfect fluid plus a 
coupling term A{Q) which characterizes a universal coupling between a fundamental scalar 
field Q and ordinary (plus dark) matter. 
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In Section 2, we have presented an explicit model for inflation, by constructing an in- 
flationary potential that, with proper choice of slope parameters, satisfies the main obser- 
vational constraints from WMAP data, including the spectral index of scalar perturbations 
and tensor-to-scalar ratio. 

In Section 3, we have first derived a set of autonomous equations, by utilizing a fun- 
damental variational principle, that in a compact form describes the evolution of different 
cosmological parameters, namely Qq, Wq, w^, e and ag, as a system of four differential 
equations, of which only three are linearly independent (cf (3.17) - (3.20)). By further 
general considerations, we have shown how the parameters q and w^^^ can be determined 
from a solution of the above system. As discussed in the body of text, the system of equa- 
tions (3.17)-(3.20) could be analytically solved only by making a reduction in the number 
of free parameters or by imposing additional constraints. In this work, one of our aims 
was to keep the model as general as possible, but for being able to find analytic solutions 
the number of parameters was restricted to four, neglecting the radiation component, and 
making a reasonable additional assumption that Q = alna + const at the present epoch. 

First by examining the case with minimal coupling, A{Q) = 1, a class of exact (ana- 
lytical) solutions has been found (cf eqs. (3.37)-(3.42)), which find interesting applications 
for the present-day cosmology. The general solution found in the minimal coupling case 
has the behavior that it is independent of the sign of a (i.e. the sign of Q). Thus the 
direction of a "rolling" scalar field Q does not seem to have any significant effect (which 
also directly followed when looking at the scalar field Lagrangian (cf eq. (3.1)), except in 
the shape of the potential. It is found that the critical value a^^^^. = 1.48 separates the 
parameter spaces of a such that a < a^^^ allows a late time acceleration while a > a^^^ 
does not. Thus the characteristic of the scalar field Q acting as an additional self-repulsive 
or self-attractive form of energy is merely determined by the magnitude of the velocity 
of the field, d{KQ) / din. a = a. In several interesting cases we have found a closed form 
expression for (reconstructed) quintessence potential V{Q). 

As the combination of WAMP and type la supernova observations show a significant 
constraint on the present-day DE equation of state, Wq = — 0.9411q5oi; for the mean 

value ujg 0.941, we require \a\ ~ 0.4207y/n^ ~ 0.36, while the WMAF-^SSS bound 

1 ^ Wq < —0.82 may be satisfied for \a\ < 0.62. Of course, a = simply represents the 
cosmological constant case {wq = — 1). Claiming the same range of — 1 < Wq < —0.82 for 
Wq at redshift z > imposes a more restrictive constraint on the slope of the potential a 
being smaller than 0.6. When looking at the evolution of different cosmological parameters 
(Jig, Q^, Wq, e, w^g, q), we find that, for smaller values of a, the model shows a late time 
accelerated expansion (for z < 1), while a matter dominance at early times. These features 
are in agreement with recent WMAP and supernova observations. 

To see how a non-minimal coupling, ^ 0, might affect the cosmic expansion, we 
studied the simplest case of an exponential coupling A{Q) oc e'^^^/'^p\ which implies 
Uq = X- III this case the solution is found to have a dependenc on the sign of the slope 
parameter a and the coupling Uq. A replacement of a by —a is found to be equivalent 
to the replacement of ag by — ctg- Moreover, a positive coupling is found to decreases the 
dark energy equation of state Wq, with respect to its value in the ag = case, while this 
effect is opposite for ctg < 0. Thus, for a fixed a, the Uq > solution could make the 
energy represented by Q more repulsive, as compared to the Og = case. The coupling 
dependence of other parameters just resemble this fact (ct^ > in our convention just 
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means a and having the same sign). For \olq\ < 0.1, and at low rcdshifts, the present- 
day values of the cosmological parameters showed almost no ag-dependence. That is, an 
observable effect on the evolution of cosmological parameters, such as w^^ and fig can be 
expected to be seen only for a strong matter-scalar coupling, like \olq\ ^ 0.1. The type la 
supernova data may favor a small value for matter-quintessence coupling, like ~ 0.06. 

We have also shown how in principle a non-minimal matter-scalar coupling can alter 
the evolution of the cosmological parameters. In general the coupling ckq always appears 
in combination with the matter density (cf eq. (3.34)). As the mass of the scalar field Q 

1 /2 

can be determined by (d^ ^ff /'^Q^) evaluated at a local minimum and the scalar-matter 
coupling in Y^^iff) can involve a p^-dependent term, the mass of a scalar field depends, 
in principle, on the ambient matter distribution. Thus in a more sophisticated model, not 
treating matter as an isotropic perfect fluid, the mass of the scalar field can vary locally 
due to a possibly strong local variation of on small scales. 
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A. Appendix: 

Corresponding to the action (3.1), the equations of motion that describe gravity, the scalar 
field Q and the background fields (matter and radiation) are given by 



These equations may be supplemented with the equation of motion of a barotropic perfect 
fluid, which is given by 




\9i.vR\-\ (V^QV.Q)+^ (VQ)'5^.+^y (Q)5m.-^^' 



(Q)r^, = 0, (A.l) 




(A.2) 




(A.3) 



Combining the (tt) and {xx) components of the equation (A.l), we get 




(A.4) 



Dividing this equation by and then using the substitution in (3.16), yields 




(A.5) 



Multiplying eq. (3.11) with Q and using the identities 



pQ = QQ^ V, Pq (1 + wq) = Q 



(A.6) 
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which follow from eq. (3.16), we get 

Pq + ^HPq (l + w^)=Q ^ E (1 - 3 P.- (A.7) 

i 

Multiplying (A. 7) by and then using equations (3.21) and (3.15)-(3.16) leads to 

2 

eq. (3.19). Further, multiplying eq. (3.14) by and then using eq. (3.21) leads to 

3^ + ^ (1 + = ^ :4(g) (1 - P, (A.8) 
Combining this equation with the identity 

n'. = ——^p'.-2en,, (A.9) 
and then using the substitutions in (3. 15)- (3. 16), finally gives equation (3.20). 
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